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Abstract Mackenzie River discharge and bathymetry effects on sea ice in the Beaufort Sea are examined 
in 2012 when Arctic sea ice extent hit a record low. Satellite-derived sea surface temperature revealed 
warmer waters closer to river mouths. By 5 July 2012, Mackenzie warm waters occupied most of an open 
water area about 316,000 km 2 . Surface temperature in a common open water area increased by 6.5°C 
between 14 June and 5 July 2012, before and after the river waters broke through a recurrent landfast ice 
barrier formed over the shallow seafloor offshore the Mackenzie Delta. In 201 2, melting by warm river waters 
was especially effective when the strong Beaufort Gyre fragmented sea ice into unconsolidated floes. The 
Mackenzie and other large rivers can transport an enormous amount of heat across immense continental 
watersheds into the Arctic Ocean, constituting a stark contrast to the Antarctic that has no such rivers to 
affect sea ice. 
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1. Introduction 

The total extent of Arctic sea ice hit a record low in the summer of 201 2 within the era of satellite observations 
spanning over three decades [Perovich et at., 2013]. A number of explanations have been suggested, 
including vulnerability to the great Arctic cyclone of August 201 2 [Simmonds and Rudeva, 201 2], sea ice 
preconditioning from winter through early summer [Devasthale et al., 2013], combined effects of the storm 
and preconditioning [ Parkinson and Comiso, 2013], and upward ocean heat transport from enhanced ocean 
mixing [Zhang et al., 2013], Here we consider a different factor — warm waters from river discharge — in 
melting sea ice in 2012. In particular, we show that the warm waters of the Mackenzie River discharged into 
the Beaufort Sea, where sea ice had been fragmented into small floes by strong winds on the Beaufort Gyre 
(BG) in early July 201 2 before the August storm occurred. Furthermore, we examine how the seafloor controls 
the recurrent pattern of sea ice and the release of the Mackenzie discharge into the Beaufort Sea. We also 
contrast the Arctic Ocean situation with that of the Southern Ocean. 


2. Sea Ice in the Beaufort Sea Study Domain in 2012 Compared to Changes 
in 1979-2011 

Sea ice has reduced severely in the Beaufort Sea [Ogi and Rigor, 201 3; Perovich et al., 2007]. The current sea 
ice cover consists of thinner, younger, and weaker ice that facilitates lead formation and influences halogen 
chemical processes [ Nghiem et al., 2013] and mercury deposition [Steffen et al., 2013]. Sea ice reduction 
in the Beaufort Sea may also shorten the snow season and adversely impact the habitat of ringed seals 
[Kelly et al., 2007], 

In this study, the Beaufort Sea domain of interest is bounded by longitudes 123°W-147°W, latitude 74°N in 
the north, and the coastline along the south side. Within this region, satellite active microwave scatterometer 
(AMS) and passive microwave radiometer (PMR) data are used to observe the change in sea ice extent (SIE) in 
early summer over the past three decades. A difficulty in summer observations of sea ice is that surface melt 
can impact PMR signatures of sea ice, causing large uncertainties in the measurement of SIE. 

To alleviate this problem, an AMS algorithm was developed to determine SIE, primarily using a quadruplet 
combination of relative AMS values rather than the absolute magnitudes of the signature and thus improving 
SIE accuracy [Nghiem et al., 2005]. QuikSCAT data in 2000-2009 are processed with the AMS algorithm to 
determine SIE average over 1-7 July in each year. Then the AMS results are applied to select an optimal value 
of PMR sea ice concentration (SIC), available at the National Snow and Ice Data Center [201 3] in Boulder, 
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Colorado, to obtain a long-term time series of 
SIE in the Beaufort Sea study domain consis- 
tently between both AMS and PMR data sets. 
SIE results in 1979-2012 for SIC above the 
optimal threshold of 0.1, which best matches 
AMS and PMR SIE in the common decade of 
2000-2009, are used to calculate open water 
extent. Figure la shows a large extent of 
open water in 2012, a larger open water 
extent in 2008, and the record high extent 
in 1998. Rawlins et at. [2009] reported a re- 
cord discharge from the Mackenzie River in 
May 1998, and Perovich et al. [201 1] showed 
significant underside ice melt resulting from 
solar heat input in 2008. Thus, we will 
examine the role of warm waters from river 
discharge in accelerating sea ice melt, 
which leads to earlier and more extensive 
open water having a lower albedo and 
allowing more solar heat absorption 
[Perovich et al., 2007]. 

The change in SIE is affected by dynamic and thermodynamic conditions of sea ice. In the Beaufort Sea, sea 
ice characteristics can change depending on the variable strength of BG activities driven by atmospheric and 
oceanic forcing. To quantify the BG strength, average values of sea ice velocity in the NW and SW orthogonal 
components at the center of the study domain in May-June are calculated from buoy data provided by the 
International Arctic Buoy Programme [Rigor, 2005], The sea ice velocity (Figure 1 b) is correlated with the open 
water extent in 1 979-2012 with a correlation coefficient of 0.77 for the NW and 0.68 for the SW component. 
When the data are detrended to assess the correlation of the interannual variability, the correlation coeffi- 
cients remain relatively high at 0.72 and 0.60 for the NW and SW components of the sea ice velocity, 
respectively. The overall ice velocity was mostly toward the west, moving away from the Canadian Arctic 
archipelago toward the Canada Basin unrestricted by land, facilitating sea ice divergence. Thus, strong BG 
activities in May-June can fragment sea ice and precondition it for enhanced melting from the subsequent 
intrusion of the Mackenzie warm waters in early July 201 2. 

3. Bathymetry Effects 

To examine further details in the case of the large open water extent in July 201 2, Moderate Resolution 
Imaging Spectroradiometer (MODIS) images are used to observe sea ice change in this region. Figure 2 is 
an overlay of two images acquired on 14 June 2012 and 17 May 2013 by the Terra MODIS. Each image is a 
composite of bands 1, 2, and 7 (620-670, 841-876, and 2105-2155 nm wavelengths) with a pixel size of 
250 m in the MODIS Rapid Response subset generated for near-real-time support of the Bromine, Ozone, 
and Mercury Experiment [Nghiem et al., 201 3]. The sea ice pattern on 14 June 2012 (light blue) reveals a 
band of landfast sea ice off the coast of the Mackenzie River Delta. The edge of this landfast sea ice traces 
consistently well along a lead (black elongated feature) opened up in the sea ice cover (dark cyan) on 
17 May 2013. 

The recurrence of the sea ice fracture suggests a control exerted by the shallow bathymetry shaping sea ice 
patterns in this littoral zone. In fact, the lead in 2013 and the landfast ice edge in 2012 are both found to 
conform the best to the 25 m isobath by a bathymetry conformity analysis [Nghiem et al., 201 2] with data 
from the International Bathymetric Chart of the Arctic Ocean (IBCAO) [Jakobsson et al., 2008, 201 2], The ice 
edges conform well to the 25 m isobath (Figure SI in the supporting information), regardless how close 
(i.e., steep slope) or how far apart (i.e., gradual slope) the other isobaths are separated from the 25 m isobath. 
These results show that sea ice patterns are controlled by bathymetric depth and are independent of seafloor 
slope (see Figure SI and Table SI). 



1980 1985 1990 1995 2000 2005 2010 

Year 


Figure 1 . (top) For open water extent within the study domain in the 
Beaufort Sea observed by satellite PMR data (1979-2012) and AMS data 
(2000-2009) averaged over 1-7 July and (bottom) for ice velocity in the NW 
and SW orthogonal components at the center of the study domain averaged 
over May-June in each year. A positive velocity value indicates ice moving 
toward the NW or the SW direction. 
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Constrained by bathymetric effects, the 
stability of the landfast ice extent has 
been observed over the past several 
decades. From Advanced Very High 
Resolution Radiometer and Landsat data 
in 1 986, Dean et at. [1 994] saw a similar 
band of landfast sea ice off the 
Mackenzie River Delta. Using Radarsat 
imagery acquired in 1996-2004 across 
the eastern Chukchi Sea to the western 
Beaufort Sea, Mahoney et al. [2007] found 
broad stable extensions of landfast sea 
ice, which did not show any significant 
difference from the 1970s. The recur- 
rence of the Mackenzie landfast ice 
feature, which determines the timing of 
the warm river discharge breakthrough 
and thus the timing of sea ice melt 
acceleration, can help in sea ice fore- 
casting for operational applications. 

It is noted that the landfast sea ice edge 
in 2012 and the lead in 2013 trace well 
along the 25 m isobath from IBCAO version 2, except over the local seafloor features marked inside the red 
ellipse (Figure SI). Those features were erroneous artifacts in IBCAO version 2, which was corrected by new 
survey data in IBCAO version 3. The new 25 m isobath from IBCAO version 3 matches all along the 2012 ice 
edge and the 2013 lead very well (Figure SI). These results suggest that recurrent sea ice patterns, once 
detected, can be useful for planning future bathymetry survey missions. 

The pattern of bathymetric isobaths can be a good indicator for landfast sea ice extent in the Beaufort Sea 
[Mahoney et al., 2007], and sea ice can conform to bathymetry recurrently in different years (Figure SI). 
However, the ice edge may not necessarily have to follow exactly along the same isobath all the time due to 
anomalous atmospheric (winds) or oceanic (currents) forcing that may become so strong to overwhelm 
bathymetry effects. 

4. Warm Waters of the Mackenzie River 

To investigate the distribution pattern of surface waters in the Beaufort Sea, two different algorithms are used 
to measure surface temperature from MODIS. The first is the ice surface temperature (1ST) algorithm (called 
AG1 hereon), which was originally developed with a split-window technique to map snow or ice surface 
temperature at a resolution of 1 km [ Hall et al., 2004]. This is the standard MODIS 1ST product, MOD29 [Riggs 
et al., 2006], implemented at the Goddard Space Flight Center and available through National Snow and Ice 
Data Center with approximately 24 h latency. Here AG1 is applied to measure temperature of the ocean water 
surface or sea surface temperature (SST). 

The second algorithm (AG2) produces daily global 1 km SST maps using Level-2 preprocessed data products 
from the Global High-Resolution SST Pilot Project and in situ ship and buoy measurements from the Global 
Ocean Data Assimilation Experiment server. A multiscale two-dimensional variational blending algorithm was 
used to combine all input data products and output SST [Chao et al., 2009]. The SST data are available at the 
Jet Propulsion Laboratory (JPL) Physical Oceanography Distributed Active Archive Center with approximately 
24 h latency [JPL OurOcean Project, 2013]. 

To observe the change caused by warm waters from the Mackenzie discharge, open water surface temper- 
ature (OWST) in the Beaufort Sea study domain is compared before and after the landfast sea ice barrier was 
breached to release the waters from the Mackenzie River into the Beaufort Sea. On 14 June 2012, while the 
open water area (Figure 3a) was estimated at about 210,600 km 2 , Mackenzie waters were still obstructed by 
the landfast sea ice barrier. This is evident in Figure 2 showing the band of landfast sea ice (light blue) across 



Blue Sea ice on both dates Cyan Sea ice on 17 May 2013 


Figure 2. Sea ice in the Beaufort Sea observed by a combination of bands 1 , 2, and 
7 from data acquired by MODIS aboard the Terra satellite. The 1 7 May 201 3 MODIS 
image is translucently overlain over the 14 June 2012 MODIS image. Light blue is 
for sea ice existing on both dates, cyan for sea ice only on 17 May 2013, and black 
for no sea ice on either date. The landfast sea ice in 2012 is observed as the light 
blue band in the littoral zone along the coast of the Mackenzie Basin, and the large 
long lead in 2013 is seen as the long black fracture tracing along the edge of the 
2012 landfast ice. 
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the mouths of streams where the warm 
waters intruded into the ice barrier and 
created several fingers of water overflows 
that had not yet reached to the landfast ice 
edge. See Figure S2 for the close-up detail 
of the water fingers; similar overflow 
patterns were observed in 1986 [Dean 
et al., 1994]. Clear sky conditions on 14 
June allow the use of MODIS satellite data 
to estimate the Beaufort Sea OWST. As 
seen in Figure 3a, OWST was low while the 
landfast ice hindered the river discharge. A 
pattern of warm waters was observed to 
the west of Banks Island, possibly from the 
inland water runoff. 

OWST on 5 July 2012 (Figure 3b) reveals 
the extrusion of warm waters in the 
Beaufort Sea when the landfast sea ice had 
disappeared and the stream mouths were 
fully opened for water discharge from the 
Mackenzie River. By that time, the area of 
open water observed by MODIS in this 
region was 316,000 km 2 , about the size of 
the country of Italy or Poland. The warmest 
waters were observed near the coast of the 
Mackenzie Delta, e.g., 1 3°C at 1 47 km, 1 0°C at 287 km, 8°C at 350 km, and 2°C as far as 456 km from the 
Mackenzie River mouth. As the warm river waters extended farther out into the Beaufort Sea, OWST became 
lower and mixed with cooler waters in a complex swirling pattern (Figure 3b). The warm waters in the west of 
Banks Island also became more extensive and started to connect to the Mackenzie warm waters mass. OWST 
was lower toward the sea ice edge in the Canada Basin where the water temperature still maintained above 
freezing temperatures. OWST patterns obtained from AG1 and AG2 are closely consistent and thus only the 
AG1 temperature is shown in Figure 3. In fact, a comparison of OWST from AG1 and AG2 in this case shows that 
the mean difference is only 0.2°C, although AG1 OWST can be lower by 0.75°C at the minimum and higher by 
1.3°C at the maximum of the temperature range compared to AG2 values. Patterns of OWST in 1998 and 2008 
reveal a similar spatial structure with warmer waters closer to the river mouth (see Figures S3 and S4). 

To obtain a quantitative assessment of the increase in OWST after the intrusion of the warm river waters into 
the Beaufort Sea, each pixel that was observed as open water in both cases on 14 June and 5 July was se- 
lected to calculate the temperature difference. The histogram of this temperature difference is presented in 
Figure 3c. The result shows a mean increase in OWST of 6.5°C between the two cases before and after the 
Mackenzie warm waters intrusion. The statistics of the OWST difference is similar to the characteristics of the 
gamma distribution, which is traditionally used to model waiting times until a termination event, such as 
death in life testing [Hogg et al., 201 2] or sea ice extinction in this case. Quantitatively, the probability density 
function of the OWST difference is determined by the gamma distribution 

f p (AT) = (AT + T c ) k ~' exp[— (AT + T c )/G]/[d k T(k)] (1) 

where T is the gamma function, ATis the OWST difference, T c = 3.0°C, k= 12.2206, and 0 = 0.740481, which are 
determined numerically by requiring that both measured and calculated values of the probability density 
have the same mean, variance, and skewness (Figure 3c). 

5. Summary and Discussions 

In summary, two different algorithms were used with MODIS satellite data to observe the intrusion of warm 
waters from the Mackenzie River discharge into the Beaufort Sea. Results from both algorithms agreed well 
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Figure 3. Open water surface temperature (OWST) in °C, measured by MODIS 
AG1 on 14 June 2012 (a) before the Mackenzie River water broke through the 
landfast sea ice barrier, and on 5 July 2012 (b) after warm waters from the 
Mackenzie discharged into the Beaufort Sea. Sea ice is indicated by the steel 
blue color on the sea surface, (c) For the measured and calculated histograms of 
the temperature difference in °C between the two cases in Figures 3a and 3b. 
The temperature difference is calculated over each pixel that is open water in 
both cases. 
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and consistently revealed the large increase in OWST after the warm river waters broke through the landfast 
sea ice barrier. The water temperature increase significantly contributed to melting sea ice that had been 
fragmented into unconsolidated ice floes by strong winds on the Beaufort Gyre. 

In the study domain, the open water extent increased by 50% in 3 weeks from 14 June to 5 July 2012, a month 
before the great Arctic cyclone occurred in August 2012. In this process, the shallow coastal bathymetry can 
exert control on the formation of the landfast sea ice barrier that obstructs the river discharge until its re- 
moval. Early model results indicated that the coastal sea ice under the riverine influence could be removed 1 
to 2 weeks in advance of areas without significant river discharge [Searcy et a!., 1 996], At larger scales, ba- 
thymetry effects on sea ice formation in the Beaufort Sea were also reported [Carmack and Chapman, 2003; 
Nghiem et at., 2012]. Since bathymetry does not change significantly over long geological time scales, it 
provides a systematic memory allowing the recurrent formation of the landfast sea ice barrier in different 
years. Such a stable geological memory from the bathymetry should be included to improve the predict- 
ability in sea ice models [Nghiem et at., 2012]. 

The Mackenzie discharge typically peaks in June and remains high in July [Woo and Thorne, 2003]. For a peak 
flow of 33,300 m 3 /s [Environment Canada, 201 3], the volume of the total discharge over the 3 week period is 
equivalent to a layer thickness of 0.19 m of warm waters across the entire open water area of 316,000 km 2 on 
5 July 2012. Note that this is only a uniform equivalence representation of the warm waters mass, while actual 
distribution of the waters can have complex mixing patterns (Figure 3). The extensive waters intrusion from 
the Mackenzie River was also identifiable by terrestrial dissolved organic matter traced by an optical proxy 
from MODIS spectral data [Fichot et at., 2013]. 

The Mackenzie River has an enormous watershed of 1.8 million km 2 with the southern extent reaching to 
52.2°N (see Figure S5 and Table S2). This watershed is primarily within the continental climate regime, and the 
heat can be intense in summer when the maximum temperature may reach 32°C around latitude 53°N (e.g., 
Edmonton, Alberta, Canada). Fresh and warm Mackenzie waters reside in the surface layer (Figure 3) with the 
attendant high thermal capacity thus contributing excessive heat to melt sea ice, most effectively when the 
sea ice cover has been fragmented by the strong BG. This was the case in the 1998 melt season when ice 
velocity was high in the BG (Figure 1 ), the Mackenzie River had a record discharge [Rawlins et at., 2009], and a 
record high of open water extent was observed in the decadal satellite data (Figure 1). 

The Mackenzie River is a highly effective conveyor that transports continental heat into the Beaufort Sea 
contributing to the overall heat budget for sea ice melting. While Mackenzie discharge had no discernable 
trend in 1 972-1 999 [Woo and Thome, 2003], an increasing trend in the mean air temperature was identified in 
station data in 1948-2007 [Turner and Gyakum, 2010], especially over the Mackenzie watershed in western 
Canada. Even assuming no temperature increase in Mackenzie River waters, the warm discharge can still 
effectively melt sea ice because Arctic sea ice, especially in the Beaufort Sea, has shifted to a new regime 
dominated by seasonal sea ice that is younger, thinner, and weaker [National Academies, 2012]. Such 
diminished sea ice becomes easily fragmented into small ice floes that increase the ice interface with warm 
waters to accelerate the melt process. 

In addition to the Mackenzie, there are a number of other large rivers that discharge into the Arctic Ocean. 
Notable are the Yukon, Ob, Yenisei, Lena, and Kolyma Rivers (see Table S2), each with its immense watershed 
under a continental climate regime providing massive discharge of warm waters into the Arctic Ocean or a 
peripheral sea to melt sea ice in spring and summer. Extensive intrusions of river runoff traced by a MODIS 
optical proxy have been reported [Fichot et ai, 2013]. A large increase from the Eurasian rivers together with a 
small decrease from North America resulted in a net discharge increase by 5.6 km 3 /yr 2 into the Arctic Ocean 
during 1964-2000 with a mean annual discharge from the mouths of 72 rivers totaled to over 2.5 trillion m 3 
by 2000 [McClelland et ai, 2006], and the 201 1 Arctic Report Card showed similar volume of total annual 
discharge from large Arctic rivers [Shiklomanov and Lammers, 2011]. This massive discharge carries an 
enormous heating power of 1.0 x 1 0 1 9 J/yr for each 1°C of the warm river waters above freezing, equivalent to 
2.5 gigaton of trinitrotoluene (TNT) per °C per year. As a comparison, Antarctic mass loss was approximately 
200 ± 1 50 Gt/yr by 2006 at a rate of — 1 4.5 ± 2 Gt/yr 2 [Rignot et ai, 2011]. When projected to 201 0, the 
Antarctic total mass loss amounts only to about 10% of the Arctic total river discharge in the annual mean. 
This is a stark contrast between the Arctic with drastic sea ice decrease and Antarctic with stable (or some- 
what increasing) sea ice since there is no such massive river discharge from the Antarctic continent, primarily 
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a frozen continent located at high latitudes and high altitudes around the South Pole. However, river dis- 
charge is not the only factor leading to the difference between Arctic and Antarctic sea ice. Other factors are 
bathymetry that guides waters distribution and circulation in the Southern Ocean, and topography on the 
Antarctic continent that controls the katabatic wind. These stable geological factors contribute to the stability 
of the Antarctic sea ice cover [Nghiem, 2013]. 

In the summer melt season, warm river waters, for which the temperatures can be higher than 10°C 
(Figure 3b), contribute directly to melting sea ice. In the fall season around the time of sea ice freezeup, 
surface waters cool while the halocline stratification insulates the surface from the deeper seawater, allowing 
more sea ice to grow. At the same time in the fall, rivers also start to freezeup, drastically reducing the river 
discharge. Thus, to be an effective insulator, the stratification needs to be persistent to maintain the surface 
layer consisting of a large mass of fresh river waters that already discharged into the Arctic Ocean earlier in 
the summer. Such maintenance of the stratification requires calm-ocean conditions without significant 
mixing throughout the summer to fall freezeup. In summer 2012, the violent storm significantly enhanced 
ocean mixing that transported ocean heat upward and further contributed to sea ice melt [Zhang eta!., 2013]. 

While this paper presents surface patterns of the Mackenzie warm waters intrusion into the Beaufort Sea, 
future hydrograph surveys using ship or buoys to measure temperature and salinity profiles across the 
Mackenzie discharge domain are necessary to determine the three-dimensional mixing for a rigorous as- 
sessment of the total heat content in warm waters available to melt sea ice. In particular, there is a severe lack 
of river water temperature profile data across river channels or river mouths where warm waters discharge 
into the Arctic Ocean. Such river temperature data are required to calculate the regional as well as total 
thermal heating power carried by the river waters to melt sea ice. This is a critical knowledge gap that needs 
to be addressed by systematic and well-calibrated networks of sensors for long-term river temperature 
measurement. While in situ data are available at river gauging stations in land, stream flow and river dis- 
charge can be measured with multidecadal data sets acquired by multiple satellite sensors [Brakenridge et at., 
2007] at many river reaches including locations near river mouths for long-term quantitative analyses. 
Rigorous interdisciplinary studies for all large Arctic rivers are necessary to understand the overall contribu- 
tion and interactions of the discharged warm waters in melting sea ice. In this respect, it is crucial to monitor 
river change, including discharge and temperature, with adequate, continuous, and accurate measurements 
by river gauges, buoy networks, and satellite sensors. 
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Introduction 

The auxiliary material includes: two tables and five figures. The first table (Table SI) is for the 
shortest distances and steepest slopes between the 25-m isobath to the 10-m isobath and to 
the 50-m isobath from IBCAO Version 3. The second table (Table S2) contains information of 
Arctic rivers. The first figure (Figure SI) shows the conformance of sea ice patterns to 
bathymetry. The second figure (Figure S2) shows fingers of open water from Mackenzie warm 
waters intruded into the band of landfast sea ice observed by MODIS. The third and fourth 
figures (Figures S3 and S4) are for surface temperature of open water in the Beaufort Sea. The 
last figure (Figure S5) presents the watershed of Mackenzie River. 

1. Table SI: Shortest distances and steepest slopes between the 25-m isobath to the 10-m 
isobath and to the 50-m isobath from IBCAO Version 3. The shortest distance is calculated as 
the radius of the smallest circle with its center located on the 25-m isobath and its perimeter 
touched the 10-m (or 50-m) isobath. The steepest slope is the ratio of the depth change 
between 25 m and 10 m (or 50 m) over the respective shortest distance. In all cases, the edge 
of landfast ice (6/14/2012) and the lead in 2013 (5/17/2013) were closely located along the 25- 
m isobath. This result illustrates that sea ice patterns conform to bathymetric depth and do not 
depend on bathymetric slope. 

1.1 Column "Longitude on 25-m isobath", degrees, west longitude on the 25-m isobath. 

1.2 Column "Distance between 25-m and 10-m isobath", km. 

1.3 Column "Slope between 25-m and 10-m isobath", m/km. 

1.4 Column "Distance between 25-m and 50-m isobath", km. 

1.5 Column "Slope between 25-m and 50-m isobath", m/km. 

2. Table S2: Watersheds of several large rivers in the Arctic. Data from ArcticRIMS, a Regional, 
Integrated Hydrological Monitoring System for the Pan-Arctic Land Mass, available on the 
Internet at http://rims.unh.edu/, accessed Dec. 2013. For comparison, the total surface area of 
Greenland is 2,166,086 km 2 . 

2.1 Column "River", river name. 

2.2 Column "Watershed", km 2 , area of watershed. 
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2.3 Column Southernmost reach", degree, north latitude of the watershed southernmost 
reach. 

2.4 Column "Sea basin", name, the sea to which the river discharge goes. 

3. Figure SI: Edge of the landfast sea ice on 14 June 2012 (white curve) and edge of the lead on 
17 May 2013 (yellow curve), observed with MODIS satellite data, overlain on the IBCAO 
bathymetry map from Version 2 (panel a) and Version 3 (panel b). The ice edges conform well 
to the 25-m isobath except over the local seafloor features marked inside the red ellipse. Those 
features were erroneous artifacts in IBCAO Version 2, and were corrected by new survey data in 
IBCAO Version 3. Sea ice patterns are thus affected by depth and independent of seafloor 
slope. 

4. Figure S2: Fingers of open water (dark blue) from Mackenzie warm waters intruded into the 
band of landfast sea ice (light blue) observed by MODIS on 14 June 2012, when the river 
discharge had not yet broken through the landfast ice barrier spanning across the Mackenzie 
Delta. 

5. Figure S3: Open-water surface temperature (OWST) in °C: Panel a for 5 July 1998, and panel 
b for 5 July 2008 (same day and month as in Figure 3 for 2012) overlain on IBCAO Version 3 
bathymetry. These results are from the Pathfinder, Level 3, Version 5, Advanced Very High 
Resolution Radiometer (AVHRR) sea surface temperature product at ftp://podaac- 
ftp.jpl.nasa.gov/OceanTemperature/avhrr/L3/pathfinder_v5/ (accessed Dec., 2013). OWST 
data are not valid in cloud-covered or sea ice areas, which appear as missing data in the maps. 

6. Figure S4: Open-water surface temperature (OWST) in °C, measured by MODIS AG1 on 5 July 
2008 (same day and month as in Figure 3 for 2012). This is similar to the AVHRR OWST map on 
the same date (auxiliary Figure S3b); however, the difference between MODIS and AVHRR cloud 
and ice masking results in some differences in missing data areas. For 5 July 1998, MODIS had 
not been launched yet and thus MODIS data were not available to obtain the 1998 OWST map. 
See auxiliary Figure S3a for OWST map from AVHRR on 5 July 1998. 

7. Figure S5: Mackenzie River watershed, represented by the blue region in North America, with 
a basin area of 1.8 million km 2 reaching south to 53°N. This map is adapted from North 
America Watersheds, Cooperative Educational Service Agency 10 (CESA 10), available on the 
Internet at http://moodle.cesal0.kl2.wi.us/, accessed Nov. 2013. 
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Table SI. Shortest distances and steepest slopes between the 25-m isobath to the 10-m 
isobath and to the 50-m isobath from IBCAO Version 3. The shortest distance is calculated as 
the radius of the smallest circle with its center located on the 25-m isobath and its perimeter 
touched the 10-m (or 50-m) isobath. The steepest slope is the ratio of the depth change 
between 25 m and 10 m (or 50 m) over the respective shortest distance. In all cases, the edge 
of landfast ice (6/14/2012) and the lead in 2013 (5/17/2013) were closely located along the 25- 
m isobath. This result illustrates that sea ice patterns conform to bathymetric depth and do not 
depend on bathymetric slope. 


Longitude 
on 25-m 
isobath 

Distance between 
25-m and 10-m 
isobath 

Slope between 
25-m and 10-m 
isobath 

Distance between 
25-m and 50-m 
isobath 

Slope between 
25-m and 50-m 
isobath 

137.16°W 

6.5 km 

2.31 m/km 

14.8 km 

1.69 m/km 

136.81°W 

19.3 km 

0.78 m/km 

15.9 km 

1.57 m/km 

134.70°W 

22.3 km 

0.67 m/km 

29.3 km 

0.85 m/km 

132.23°W 

18.7 km 

0.80 m/km 

51.4 km 

0.49 m/km 

131.37°W 

16.1 km 

0.93 m/km 

54.7 km 

0.46 m/km 

129.38°W 

34.8 km 

0.43 m/km 

40.7 km 

0.61 m/km 


Table S2. Watersheds of several large rivers in the Arctic. Data from ArcticRIMS, a Regional, 
Integrated Hydrological Monitoring System for the Pan-Arctic Land Mass, available on the 
Internet at http://rims.unh.edu/, accessed Dec. 2013. For comparison, the total surface area of 
Greenland is 2,166,086 km 2 . 


River 

Watershed 

Southernmost reach 

Sea basin 

Mackenzie 

1,783,912 km 2 

52.2°N 

Beaufort Sea 

Yukon 

833,232 km 2 

58.8°N 

Bering Strait 

Ob 

2,994,238 km 2 

45.3°N 

Kara Sea 

Yenisei 

2,537,404 km 2 

45.7°N 

Kara Sea 

Lena 

2,460,742 km 2 

52.2°N 

Laptev Sea 

Kolyma 

651,631 km 2 

60.6°N 

East Siberian Sea 


S3 


Supporting Information, Geophysical Research Letter, 2014 



Figure SI. Edge of the landfast sea ice on 14 June 2012 (white curve) and edge of the lead on 17 
May 2013 (yellow curve), observed with MODIS satellite data, overlain on the IBCAO 
bathymetry map from Version 2 (panel a) and Version 3 (panel b). The ice edges conform well 
to the 25-m isobath except over the local seafloor features marked inside the red ellipse. Those 
features were erroneous artifacts in IBCAO Version 2, and were corrected by new survey data in 
IBCAO Version 3. Sea ice patterns are thus affected by depth and independent of seafloor 
slope. 
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Figure S2. Fingers of open water (dark blue) from Mackenzie warm waters intruded into the 
band of landfast sea ice (light blue) observed by MODIS on 14 June 2012, when the river 
discharge had not yet broken through the landfast ice barrier spanning across the Mackenzie 
Delta. 
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Figure S3. Open-water surface temperature (OWST) in oC: Panel a for 5 July 1998, and panel b 
for 5 July 2008 (same day and month as in Figure 3 for 2012) overlain on IBCAO Version 3 
bathymetry. These results are from the Pathfinder, Level 3, Version 5, Advanced Very High 
Resolution Radiometer (AVHRR) sea surface temperature product at ftp://podaac- 
ftp.jpl.nasa.gov/OceanTemperature/avhrr/L3/pathfinder_v5/ (accessed Dec., 2013). OWST 
data are not valid in cloud-covered or sea ice areas, which appear as missing data in the maps. 
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Figure S4. Open-water surface temperature (OWST) in oC, measured by MODIS AG1 on 5 July 
2008 (same day and month as in Figure 3 for 2012). This is similar to the AVHRR OWST map on 
the same date (auxiliary Figure S3b); however, the difference between MODIS and AVHRR cloud 
and ice masking results in some differences in missing data areas. For 5 July 1998, MODIS had 
not been launched yet and thus MODIS data were not available to obtain the 1998 OWST map. 
See auxiliary Figure S3a for OWST map from AVHRR on 5 July 1998. 
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Figure S5. Mackenzie River watershed, represented by the blue region in North America, with a 
basin area of 1.8 million km2 reaching south to 53oN. This map is adapted from North America 
Watersheds, Cooperative Educational Service Agency 10 (CESA 10), available on the Internet at 
http://moodle.cesal0.kl2.wi.us/, accessed Nov. 2013. 
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